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Abstract 

We present the elliptic flow of inclusive J /if/ measured in the p + p~ channel at forward rapidity 
(2.5 < y < 4.0), down to zero transverse momentum, in Pb-Pb collisions at -Js m = 2.76 TeV 
with the ALICE muon spectrometer. The pj dependence of J/if/ v'2 in non-central (20%-60%) 
Pb-Pb collisions at -\/s NN = 2.76 TeV is compared with existing measurements at RHIC and 
theoretical calculations. The centrality dependence of the px-integrated elliptic flow, as well as 
the pj dependence in several finer centrality classes is presented. 



1. Introduction 

Charmonium production in heavy ion collisions has been studied at different energies and 
with different collision systems, ever since the J/if/ suppression induced by color screening of its 
constituent quarks was proposed as a signature of the formation of a quark gluon plasma (QGP) 
in heavy-ion collisions [1]. The recent measurement of the J/if/ production in Pb-Pb collisions at 
forward rapidity performed by ALICE at the LHC |2] clearly showed less suppression compared 
with SPS and RHIC results 0,13]. At RHIC energies, the preliminary result from the STAR 
collaboration showed a 1/if/ elliptic flow in Au-Au collisions at ys NN = 200 GeV J5[] consistent 
with zero within uncertainties in the measured pj range (0-10 GeV/c). The measurement of 
quarkonium elliptic flow is especially promising at the Large Hadron Collider (LHC) where the 
high energy density of the medium and the large number of cc pairs produced in Pb-Pb collisions 
is expected to favor the flow development and regeneration scenarios. 

2. Data analysis and results 

The ALICE detector is described in 16J. At forward rapidity (2.5 < y < 4.0) the production 
of quarkonium states is measured in the muon spectrometer down to pj = 0. The data sample 
used for this analysis corresponds to 17 M dimuon unlike sign (MU) triggered Pb-Pb collisions 
collected in 201 1 . It corresponds to an integrated luminosity L; nt * 70 jub -1 . The event and muon 
track selection are the same as described in [7], except for an additional requirement of the event 
vertex position |Z vtx | < 10 cm to ensure a flat event plane distribution. J/tf/ candidates are formed 
by combining pairs of opposite-sign (OS) tracks reconstructed in the geometrical acceptance of 
the muon spectrometer. 
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Figure 1 : Left: V2 extraction with the event plane method in which the J/i// raw yield is plotted as a function of Atp using 
a fit to the data with dN'^/dAip = A X (1 + 2v° bs cos 2A(p). Right: Inclusive ]/ifi v 2 in the centrality bin 20%-60% as 
a function of p~[. A comparison to STAR results and to two parton transport model calculations IIOL fllll are shown. 
The vertical bars show the statistical uncertainties, and the boxes indicate the point-to-point uncorrelated systematic 
uncertainties, which are dominated by the signal extraction. 



The J /iff V2 is measured using event plane based methods |8]. To make a direct comparison 
with lower energy measurements, the inclusive J/if/ V2 (pj) was calculated in the same centrality 
range 20%-60% as at RHIC [5], as discussed in detail in The V2 is extracted by fitting the J /iff 
raw yield as a function of Aip = ^dimuon - ¥ep,2 with dN ] ^/dAip = A x (1 + 2v'2 bs cos 2A<p), where 
A is a normalization constant (the standard event plane method), as shown in Fig. [TJ(left panel). 
At LHC, the event-plane-resolution-corrected v-i of J /iff with 2 < pj < 4 GeV/c is different from 
the STAR preliminary measurement which is compatible with zero in all the measured pj range, 
as shown in Fig. Q] (right panel). Two model calculations based on transport mechanism which 
include a J /if/ regeneration component from deconfined charm quarks in the medium iTToL [Till 
are compared with data. These two models differ mostly in the rate equation controlling the J/if/ 
dissociation and regeneration. In both models about 50% of the produced J /if/ mesons originate 
from regeneration in QGP in the most central collisions. On one hand, thermalized charm quarks 
in the medium will transfer a significant elliptic flow to regenerated J/if/. The maximum \>2 at 
pj ~ 2.5 GeV/c results from a dominant contribution of regeneration at lower pj with respect to 
the initial J /if/ component. Both models are able to qualitatively describe J/if/ V2(Pt) data as both 
were also able to describe the earlier R\\ measurement 12J]. 

Consistent results in the same centrality bin are obtained with an invariant mass fit technique, 
in which we fit the V2 = (cos2Ai^) vs. invariant mass (m w ) as described in 01211 . The method 
involves calculating the v-i of the OS dimuons as a function of m w and then fitting the resulting 
V2 (m-fj/j) distribution using: viim^) = v^aim^) + v^im^ll - a(m m )], where v" 8 is the J/iff 
elliptic flow and v 2 6 is the background flow (parametrized using a second order polynomial 
function in this analysis). a(m m ) = S/(S + E) is the ratio of the signal over the sum of the 
signal plus background of the m ff distributions, aim^) is extracted from fits to the OS invariant 
mass distribution in each pj and centrality class. The OS dimuon invariant mass distribution 
was fitted with a Crystal Ball (CB) function to reproduce the J /if/ line shape, and either a third 
order polynomial or a variable width gaussian to describe the underlying continuum. The CB 



function connects a Gaussian core with a power-law tail II 1 3fl at low mass to account for energy 
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Figure 2: Left: (pr) 11 "™ 1 and v% extraction with the fit invariant mass technique. Right: Event plane resolution corrected 
J/ij/ V2 as a function of centrality of ]/if/ with pj > 1.5 GeV/c. The vertical bars show the statistical uncertainties, and the 
boxes indicate the point-to-point uncorrected systematic uncertainties, which are dominated by the signal extraction. 



49 loss fluctuations and radiative decays. The combination of several CB and underlying continuum 

so parametrization were tested to assess the signal and the related systematic uncertainties. The J/if/ 

51 V2 in each pj and centrality class was determined as the average of the v" g obtained by fitting 

52 V2(w wu ) with various background shapes, while the corresponding systematic uncertainties were 

53 defined as the r.m.s. of these results. A similar method is used to extract the uncorrected average 

54 transverse momentum (/?T) uncor of the reconstructed J/if/ in each centrality and pj class. Fig. [2] 

55 (left panel) shows typical fits of the OS invariant mass distribution (top left), the (cos 2(0— *Pep,2)) 

56 (bottom left) and (pT) uncor (middle left) as a function of in the 20%-60% centrality class. The 

57 obtained J /if/ (/?t) ui1C01 is used to locate the ALICE points when plotted as a function of transverse 

58 momentum. 

59 Fig. [2] (top right) shows V2 for inclusive J/i// with pj > 1.5 GeV/c as a function of centrality. 
eo The vertical bars show the statistical uncertainties while the boxes indicate the point-to-point 

61 uncorrected systematic uncertainties from the signal extraction. The measured V2 depends on 

62 the pi distribution of the reconstructed i/if/. Therefore, (/?T) uncor of the reconstructed V2 is also 

63 shown in Fig. [2] (bottom right) as a function of centrality. For the two most central bins, 5%-20% 

64 and 20%-40% the inclusive J /if/ v 2 for p T > 1.5 GeV/c are 0.101 ± 0.044(stat.) + 0.003(syst.) 

65 and 0.116 + 0.045(stat.) ± 0.017(syst.), respectively. For the two most peripheral bins the V2 is 
ee consistent with zero within uncertainties. Although there is a small variation with centrality, the 
67 (pr) lmcor stays in the range (3.0, 3.3) GeV/c indicating that the bulk of the reconstructed J /if/ are 
ea in the same intermediate pj range for all centralities. Thus, the observed centrality dependence 

69 of the V2 for inclusive J/if/ with pj > 1 .5 Ge V/c does not result from any bias in the sampled pj 

70 distributions. 

71 Fig.|3]shows the inclusive J /if/ V2(pr), using the invariant mass fit technique, for central, semi- 

72 central and peripheral Pb-Pb collisions at 2.76 TeV. In the semi-central (20%-40%) case, taking 

73 into account statistical and systematic uncertainties, the combined significance of a non-zero v% 

74 in 2 < pi < 6 GeV/c range is 3cr. At lower and higher transverse momentum the inclusive 

75 J/0- V2 is compatible with zero within uncertainties. In most central (5%-20%) and peripheral 

76 (40%-60%) case, the large uncertainties do not allow any firm conclusion. 
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Figure 3: J/i/< V2 as a function of px in various centrality bins: 5%-20%, 20%-40% and 40%-60%. The vertical bars 
show the statistical uncertainties, and the boxes indicate the point-to-point uncorrelated systematic uncertainties from the 
signal extraction. 

3. Summary and conclusion 

In summary, we reported the ALICE measurement of inclusive J/i/r V2 at forward rapidity in 
Pb-Pb collisions at ys NN = 2.76 TeV. For non-central (20%-60%) collisions a hint of a non-zero 
J/if/ elliptic flow is observed in the intermediate pj range in contrast to the zero V2 observed at 
RHIC. Indication of a non-zero ]/iff is also observed in semi-central (20%-40%) collisions at 
intermediate pj. The integrated V2 of }/tfr with pi > 1.5 GeV/c in 5%-40% collisions also shows 
a non-zero behavior. These measurements complement our earlier results on J/i/r suppression, 
where a smaller suppression was seen at low transverse momentum at the LHC compared to 
RHIC |0, 0. [3. Both results taken together could indicate that a significant fraction of the 
observed J/tfr are produced from a (re)combination of the initially produced charm quarks. Our 
J/i/r elliptic flow results in Pb-Pb collisions at ys NN = 2.76 TeV are in qualitative agreement with 
transport models that are able to reproduce our i/ift /?aa measurement. 
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